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stable equilibrium configurations. Considering the transparency condition, we study 
the stability of thin shells around black holes, showing that our analysis is in agreement 
with previous results. Applying the analysis to traversable wormhole geometries, by 
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that the latter may be significantly increased by considering appropriate choices for 
the redshift function. 
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1. Introduction 

The study of hypersurfaces of discontinuity plays a fundamental role in general relativity. 
Pioneering work can be traced back to Sen [1 , Lanczos Darmois [3], and later 
by O'Brien and Synge [Hj, Lichnerowicz t 5j, and Israel 0|, amongst others. Several 
approaches to deduce the geometric conditions of the surface layer were carried out. 
For instance, Sen in a relatively unknown paper [JJ, used continuous four- dimensional 
coordinates across the layer, an approach later applied by Lichnerowicz, and O'Brien 
and Synge. Darmois expressed the surface properties of the discontinuity of the extrinsic 
curvature across the layer as a function of the surfaces intrinsic coordinates, which was 
later generalized by Israel jSj, using independently defined four-dimensional coordinate 
systems x± of the two different manifolds glued at the junction surface. Since then, the 
hypersurfaces of discontinuity have had an extensive range of applicability. It is perhaps 
important to emphasize that the matching of an interior dust solution to an exterior 
Schwarzschild spacetime, by Oppenheimer and Snyder, provided the first insights of the 
gravitational collapse into a black hole [7j. One can mention more recent applications of 
the thin shell formalism, namely, gravitational collapse of radiating shells, the evolution 
of bubbles and domain walls in cosmological settings and in inflationary models, the 
structure and dynamics of voids in the large scale structure of the universe, shells around 
black holes and their respective stability, signature changes, matchings of cosmological 
solutions, and applications to the Randall-Sundrum brane world scenario to wormhole 
physics. An analysis in thick gravitating walls expanded in powers of the thickness of 
the wall has also been carried out [HI Ej- The lightlike limit was analyzed in Ref. "TO], and 
an equivalence between the Darmois-Israel formalism and the distributional method for 
thin shells was established in Ref. |[ITj. Due to the extensive applications of the Darmois- 
Israel formalism, a computer algebra system was implemented to aid relativists in the 
evaluation of junction conditions and the parameters associated with thin shells "T2"j . 

The application of the Darmois-Israel formalism to black holes and wormhole 
physics is of a particular interest. Relatively to black holes, in an interesting paper 
[T""] . the Schwarzschild metric inside the event horizon was matched to the de Sitter 
solution at a spacelike junction surface. The respective properties of the layer were 
explored, and it was found that instead of a singularity, a closed world can be formed 
inside the black hole. In Ref. [T"I], a thin shell was constructed around a black hole 
and its characteristics were explored by imposing the energy conditions, and in Ref. 
[To] the respective stability of the thin shell was analyzed against spherically symmetric 
perturbations about a static solution. Relatively to wormhole physics, using the cut- 
and-paste technique, Visser constructed thin-shell wormholes, and considered a partial 
stability analysis by imposing specific equations of state [Ji"J[T^[Tj3[T""!]. In Ref. [211] thin- 
shell Schwarzschild wormholes were considered and the respective stability to spherically 
symmetric perturbations around static solutions was analyzed, in the spirit of ["""""]• It 
was later found that the inclusion of a charge [21] and of a positive cosmological constant 
[2*2] significantly increases the stability regions. Specific wormhole solutions were also 
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constructed by matching an interior traversable wormhole solution to exterior vacuum 
spacetimes at a junction interface [2H1 1231 123 I2H] ■ Thin shells with a zero surface energy 
density were analyzed in Ref. |2Sj, dust shells in Ref. [23, generic surface stresses in 
Ref . [22] , and a similar analysis for plane symmetric traversable wormholes in an anti-de 
Sitter background was extensively studied in Ref. [2E]. 

Ishak and Lake developed a formalism to analyze the stability of spherically 
symmetric thin shells by imposing spacetimes that satisfy the transparency condition 
[2Zj, which amounts to considering solutions that do not contribute with the 
discontinuity flux term in the conservation identity. The dynamics of timelike spherical 
thin shells, satisfying the transparency condition were also analyzed in Ref. [28J. In this 
work, we include this momentum flux term, which severely complicates the analysis, 
and following the Ishak-Lake approach, we analyze the stability of generic spherically 
symmetric thin shells to linearized perturbations about static equilibrium solutions. 
We then deduce a master equation dictating the stability regions. Considering the 
transparency condition, we study the stability of thin shells around black holes, showing 
that our analysis is in agreement with previous results. We also apply the general 
formalism to specific wormhole solutions, and by considering particular form functions, 
we verify that the stability regions can be significantly increased by taking into account 
appropriate choices of the redshift function. In this context, a specific application of 
the general formalism developed in this paper has been applied to traversable wormhole 
solutions [2H] supported by phantom energy [3Tj[ I3T]. 

The plan of this paper is as follows: In Section II, we outline a brief summary of 
the Darmois-Israel formalism, which shall be used throughout the work. In Section III, 
we present two generic spherically symmetric spacetimes matched together at a junction 
interface and deduce the surface stresses of the thin shell. Based on the inclusion of 
the momentum flux term in the conservation identity, fundamental relationships for 
analyzing the stability regions are also deduced. From the equation of motion, a master 
equation, dictating the stable equilibrium configurations, is deduced. In Section IV, we 
shall apply the formalism developed to thin shells around black holes. In Section V, 
the formalism is applied to traversable wormholes, imposing several choices for the form 
and redshift functions. Finally, in Section VI, we conclude. 

2. The Darmois-Israel formalism 

Consider two distinct spacetime manifolds, M.+ and Ai-, with metrics given by g~^ u {x+) 
and guu( x -)i i n terms of independently defined coordinate systems x+ and x^L. The 
manifolds are bounded by hypersurfaces S + and £_, respectively, with induced metrics 
gfj and gZ. The hypersurfaces are isometric, i.e., = = 1X1 t erms of the 

intrinsic coordinates, invariant under the isometry. A single manifold M. is obtained by 
gluing together Jvt + and M._ at their boundaries, i.e., M. = with the natural 

identification of the boundaries £ = S + = £_. In particular, assuming the continuity of 
the four-dimensional coordinates x± across S, then g~ = g+ is required, which together 
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with the continuous derivatives of the metric components dg^ v /dx a \^ = dg^ u /dx a \ + , 
provide the Lichnerowicz conditions [3]. 

The three holonomic basis vectors eu\ = d/d^ 1 tangent to E have the following 
components e£J± = <9x±/<9£*, which provide the induced metric on the junction surface 
by the following scalar product 

9ij = e (i) ■ e U) = 9^ e (i) e U)\±- (!) 

We shall consider a timelike junction surface E, defined by the parametric equation 
of the form = 0. The unit normal 4— vector, n M , to E is defined as 

df df 



± 



9 



dx a dxP 



1/2 Q f 

(2) 



dx^ ' 



with n^n^ = +1 and n^e^ = 0. The Israel formalism requires that the normals point 
from M- to M+ [Ej. 

The extrinsic curvature, or the second fundamental form, is defined as = 

r W e (i) e (i)' or 

(O^x^ 1 i 0x a dx^ \ 

Note that for the case of a thin shell Kij is not continuous across E, so that for 
notational convenience, the discontinuity in the second fundamental form is defined as 
= — K~j. In particular, the condition that = g^, together with the continuity 
of the extrinsic curvatures across E, = K^, provide the Darmois conditions [3]. 

Now, the Lanczos equations follow from the Einstein equations for the hypersurface, 
and are given by 

s) = ~ (**, - &A) > ( 4 ) 

where 5 l - is the surface stress-energy tensor on E. 

The first contracted Gauss-Kodazzi equation or the "Hamiltonian" constraint 

G^nV = \ (K 2 - KijK ij - 3 R) , (5) 

with the Einstein equations provide the evolution identity 

S ij Kij = - [T^n»} + _ . (6) 

The convention [X]± = X+| E - X"| s and X = (X+| E + X~\ s )/2 is used. 

The second contracted Gauss-Kodazzi equation or the "ADM" constraint 

G^n" = K\. - K n , (7) 

with the Lanczos equations gives the conservation identity 

(8) 

In particular, considering spherical symmetry considerable simplifications occur, 
namely K % j = diag (k t t , K e e , n e ^j . The surface stress-energy tensor may be written 



>=>j\i 
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in terms of the surface energy density, a, and the surface pressure, V, as = 
diag(— cr, V, V). The Lanczos equations then reduce to 

- = -^4, 0) 

V = ^{k\ + k%). (10) 

3. Generic dynamic spherically symmetric thin shells 

3.1. Junction conditions 

We shall consider, in particular, the matching of two static and spherically symmetric 
spacetimes given by the following line elements 

ds 2 ± = - e 2a±(r±) dt\ + e 2/3±(r±) dr 2 ± + r 2 ± (d6 2 ± + sin 2 6± #|) , (11) 

of M.±, respectively Using the Einstein field equation, = 8ttT^, in an orthonormal 
reference frame, (with c — G — 1) the stress-energy tensor components are given by 

P (r) = - — {2(5' r + e 2 ^ - l) , (12) 

1 p-2/3 

Pr (r) = -—(2a>r-e^ + l), (13) 

1 e~ 2<i r i 

p t (r)=- -f3' + a' + ra" + r(a') 2 -rap' , (14) 

87r r L - 1 

where we have dropped the ± subscripts as not to overload the notation, and in which 

p{r) is the energy density, p r (r) is the radial pressure, and Pt(r) is the lateral pressure 

measured in the orthogonal direction to the radial direction. 

The energy conditions will play an important role in the analysis that follows, so 

we will at this stage define the null energy condition (NEC). The latter is satisfied if 

T^ u k u > 0, where T^ v is the stress-energy tensor and any null vector. Along the 

radial direction, with k^ = (1, ±1,0,0) in the orthonormal frame, we then have the 

following condition 

F k° = p{r) + p r (r) = ^- e' 2/3 (a' + /?')> . (15) 

47rr 

The single manifold, A4, is obtained by gluing M. + and M._ at E, i.e., at 
f(r, t) — r — a(r) = 0. In order for these line elements to be continuous across the 
junction, we impose the following coordinate transformations 

e a-{a) g/3-(a) 

The intrinsic metric to E is thus provided by 

ds% = -dr 2 + a(r) 2 (d6 2 + sin 2 9 d<p 2 ) . (17) 

The imposition of spherical symmetry is sufficient to conclude that there is no 
gravitational radiation, independently of the behavior of the junction surface. The 
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position of the junction surface is given by x^(t,9, 
respective 4-velocity is 



U£ = (V ±(a) - a±(a) ^e-^W + d 2 , d, 0, 



(t(r), a(r), 9, 0), and the 

(18) 



where the overdot denotes a derivative with respect to r. 

The unit normal to the junction surface may be determined by equation (J2J) or by 



the contractions, U^n^ = and n^n^ 



-1, and is given by 



3 /8±(a)-a±(a) 



a, Ve-^W + d 2 ,0,0 



(19) 



Using equation (J3J), the non-trivial components of the extrinsic curvature are given 



by 



+ d 2 , 



K 



T ± 



a'± (V 2/3 ± + d 2 ) + a + (3' ± a 2 
y/e-W* +a? 



(20) 



(21) 



The Einstein equations, eqs. ©-(HUJ), with the extrinsic curvatures, eqs. (}2T)j) - (j2~T]) . then 
provide us with the following expressions for the surface stresses 

1 



a = — 



Ana 



e-W + a 2 



V 



■ina 



(1 + aa') (V 2/3 + d 2 ) + ad + (3 1 a a 
Ve- 2 ? + a 2 



(22) 



(23) 



If the surface stress-energy terms are zero, the junction is denoted as a boundary surface. 
If surface stress terms are present, the junction is called a thin shell. Note that the 
surface mass of the thin shell is given by m s = 4na 2 a 



3.2. Conservation identity 

Taking into account the transparency condition, [G^U^ n u ]_ = 0, the conservation 
identity, equation (jHJ), provides the simple relationship 



or 



a = -2-(a + V) 
a 



diaA) m dA 
dr dr 



(24) 



(25) 



where A = Ana 2 is the area of the thin shell. The first term represents the variation 
of the internal energy of the shell, and the second term is the work done by the shell's 
internal force. 

Now, taking into account the momentum flux term in equation ((? 
[T^U" n"]±, we have 



(Tii + T^e^dx/e- 2 ? + a 2 



(26) 
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where and Tf? are the stress-energy tensor components given in an orthonormal basis. 
The flux term, equation ([26)1 . corresponds to the net discontinuity in the momentum 
flux = U u which impinges on the shell. Thus in the general case, the conservation 
identity provides the following relationship 



a' = (a + V) + E , 

a 

where for notational convenience, we have defined H as 

1 



2ii 



e- 2 ? + a? 



Ana 



(«' + f3') ^Je-W + a? 



(27) 



(2f 



Note that this flux term vanishes in the particular case of p = —p. 
Taking into account the following relationship 



a + V 



1 



■$na 



(aa' — 1) e 2/3 + a(a' + (3')d 2 — a 2 + ad 



Ve^ + a- 

and the definition of 5, equation (J27)) takes the form 

;i + af3') e~ w + d 2 - ad 



a 



Ana 2 



Ve-W + a- 



(29) 



(30) 



which can also be obtained by taking the radial derivative of equation (|22|). Evaluated 
at the static solution ao, with a = a = 0, this reduces to 

„ , 1 



Anal 



e-^l + aoP') 



(31) 



Now, using m s = 47ra 2 cr, and taking into account the radial derivative of a', equation 
f!27jl can be rearranged to provide the following relationship 



m 



_ 1 =T-4 7 ra / r?, 
2a J '' 

with the parameter 77 defined as 77 = V'/cx', and T given by 
T = — (a + V) + 2vra S' . 



(32) 



(33) 



The physical interpretation of 77 is extensively discussed in |2*m 12*7] , and ^/rj is normally 
interpreted as the speed of sound. Equation (J32)) will play a fundamental role in 
determining the stability regions of the respective solutions which will be analyzed 
further ahead. 

For self-completeness, we shall also add the expression of H', given by 



Ana 2 



[a{a" + 13") - {a' + (3')] V e~ 2 ? + a 2 + 



a(a' + (3')(-(3'e- 2 P + a) 



Ve~ Z2 ~f f T 



a 2 



which evaluated at the static solution, ao, reduces to 

1 



"0 



Anal 



[a (a" + 0") - (a' + (3')(1 + a f3')} e~ p 



(34) 



(35) 
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3. 3. Equation of motion 

Rearranging equation f!22[) into the form 



e" 2 ^ + a 2 



e~ 2/3 - + a 2 — 47ra a , 
we deduce the thin shell's equation of motion, i.e., 

a 2 + V(a) = . 
The potential V(a) is given by 



V{a) 



F(a) 



m s \ 2 /aG\ 2 



2a J 



(36) 
(37) 

(38) 



\m s J 

where m s (a) = Ana 2 a is the mass of the thin shell. The factors F(a) and G(a 
introduced for computational convenience, are defined by 

\ ( e " 2 "- + e 
+ 1 



F(a) = e- 2 ? 

1 



G(a) = - 



2 

-2/3' 



(< 



-2/3+ 
-2/?_ 



— e 



-2/3 H 



(39) 
(40) 



respectively. 

Linearizing around a static solution situated at ao, we consider a Taylor expansion 
of V(a) around a to second order, given by 



V(a) = V(a ) + V'(a Q )(a - a ) + -V"(a )(a 



a ) 2 + O[(a-a ) 3 }. (41) 



Note that one presumes that V(a ) = 0, which is imposed in order to make the expansion 
consistent. The first and second derivatives of V(a) are, respectively, given by 



V'la) — F' — 2 



V"(a) = F" 



m s 
2a 



2a 



m s 
2a 



2 ^ ^ -2 



m s / \m s J 
aG" 



2a J \2a 



2 ( — ^ 



(42) 



(43) 



Evaluated at the static solution, at a = a , through a long and tedious calculation, 
we find V(a ) = and V'(ao) = 0. Thus, the potential, equation (}4~Tj) . reduces to 
V(a) = \V"(ao)(a — a ) 2 + 0[(a — ao) 3 ]. From the condition V'(a Q ) = 0, one extracts 
the following useful equilibrium relationship 



frn s y 
V2a 



r= [il 



F' — 2 



m. 



(44) 



Note that if V'(ao) < is verified, the potential V(a) has a local maximum at ao, 
where a small perturbation in the surface radius will produce an irreversible contraction 
or expansion of the shell. Therefore, the solution is stable if and only if V(a) has a local 
minimum at a and V"(a ) > is verified. The latter stability condition takes the form 



2a 

where \1/ is defined as 
F" 
~2 



2a 



aG 

m s 



aG 



aG" " 

m s . 



(45) 



(46) 
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3.4- The master equation 

Substituting equation (J32|) into equation (}4"5|) . one deduces the master equation given 
by 

a' m s 770 > 6 , (47) 
where r/o = r](a ) and 6 is defined as 

Now, from the master equation we deduce that the stable equilibrium regions are 
dictated by the following inequalities 

?7o>e, if a'm s >0, (49) 
770 < 6, if a'm s < , (50) 

with the definition 

5=-^-. (51) 
o~ m s 

In the specific cases that follow, the explicit form of G may become extremely messy, 
so that as in |27 , we find it more instructive to show the stability regions graphically. 

4. A dynamic shell around a black hole 

Consider that the interior and exterior spacetimes are given by equation (fTTj) . with the 
following definitions 

e 2a ± = e -2/3 ± = L_ (52) 

with m + = M and m_ = m. These are separated by a thin shell with the surface 
stresses given by 



° = -4^\f-- + a2 -f-- + a2 )- (53) 

1 / 1 - ^ + a 2 + ad l-^ + fl 2 + oa , 

V = — — — a — (54) 

87ra I Ji-m + a? Ji-^ + a 2 ' 



with the junction radius at a > 2M. 

For this case the transparency condition holds, as 2$ = T^f = and consequently 
H = 0, so that the conservation identity provides the following simple relationship 

a' = —(<t + V). (55) 
a 

Taking into account equations (|5^) - (|54*jl . then equation (|53j) takes the form 
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and at the static solution a reduces to 

i / 1 — — 1 — — \ 

a'(a ) = * , a ° - , °° , (57) 

which plays a fundamental role in determining the stability regions. 
Evaluated at the static solution, a , equations (|53 )l -(|53 j) reduce to 

<>= - 7^r(\'- — -» i- — )• (•> 




73 = a ° - a » . (59) 

87ra 1 " w ! ' : ' 

The weak energy condition (WEC) holds if a > and a + V > are satisfied, and 
by continuity implies the null energy condition (NEC), a + V > 0. The strong energy 
condition (SEC) is satisfied if a + V > and o + 2V > 0. Note that if m < M, then the 
surface energy density is positive, o > 0, and it can also be shown that the tangential 
surface pressure is also positive, V > 0. Thus, if m < M, the NEC, WEC and the SEC 
are readily satisfied. If m > M, then the surface energy density and the surface pressure 
are negative, consequently violating the energy conditions. 

The equation of motion is of the form a 2 + V(a) = 0. Using the factors of equations 
fl39>(EE0|) . we have 

m + M „ M - m , s 

F = 1 — and G = , (60) 

a a 

and the potential is given by 

M + m m 2 s {M-mf 
V(fl) = 1 -^— (61) 
If m < M, then m s (a) > and a'(a) < 0. If m > M, then m s (a) < and o~'(a) > 0. 
For both cases the master equation o'm s r] > 9 is satisfied, so that the stability regions 
supersede the energy conditions. Thus, equation pg|). dictates the following regions of 
stability 

^o<e, (62) 
which is shown below the surface plotted in figure ^ This is in agreement with previous 
results PUSH- 

5. A dynamic shell around a traversable wormhole 

5.1. Interior and exterior solutions 

Consider that the interior and exterior spacetimes are given by equation with the 
following definitions 

e 2a + = e -2 P+ = f 1 _?M\ (g3) 



r-4 



a. 



$(r_) , e- 2 ^ = (l- ^] . (64) 
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Figure 1. The regions of stability are situated below the surface, which is given by 
0. This is in agreement with previous results. See the text for details. 



We have considered that the exterior spacetime is given by the Schwarzschild solution. 
We shall once again drop the (±) subscripts, as not to overload the notation. The 
spacetime with the definition (|64j) describes a wormhole geometry f3E2], where $(r) 
and b(r) are arbitrary functions of the radial coordinate, r, denoted as the redshift 
function and the form function, respectively j22] • Note that both solutions are matched 
at a junction interface, r = a, situated outside the event horizon, i.e., a > r b = 2M, to 
avoid a black hole solution. 

Relatively to the interior solution, the field equations, eqs. fll2|) - (fTij) . with the 
definition (JMj) provide the following stress-energy scenario 

(65) 



p(r) 



57T r z 



Pr\T) = 7T 



Pt[r) 



1 


b 












-] 


8^ 






r 




A 2 



b'r — b 



b'r — b 



+ 



(66) 
(67) 



2r(r — b) 2r 2 (r — b) 

in which p(r) is the energy density; p T {r) is the radial pressure; and p(r) is the transverse 
pressure. 

Wormhole spacetimes necessarily violate the null energy condition (NEC) at the 
throat. The NEC as defined by equation ()15p now takes the form 



1 ^ylV h, 



p(r) +p r (r) 



1 



57T 



b'r -b 



2 1 




(6* 



Taking into account the flaring out condition of the throat deduced from the 
mathematics of embedding, i.e., (6 — b'r) /2b 2 > 0, evaluated at the throat b(r ) = r = r , 
and due to the finiteness of $(r), we verify that equation (jHEJ) is necessarily negative, 
i.e., Tj 1 ( / k' 1 k v < 0. Matter that violates the NEC is denoted as exotic matter. 

Recently, Visser, Kar and Dadhich, by implementing the notion of the "volume 
integral quantifier", showed that the interior wormhole solution may be supported by 
arbitrarily small quantities of averaged null energy condition (ANEC) violating matter 
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[3*3*] . although the NEC and WEC are always violated for wormhole spacetimes. In 
the spirit of minimizing the usage of the exotic matter, regions where the surface 
stress-energy tensor obeys the energy conditions at the junction were found """"3 E3J 
Recently, it was also shown that traversable wormholes may be supported by phantom 
energy "3~"""[ l3*T] , a null energy condition violating cosmic fluid responsible for the present 
accelerated expansion of the universe. In Ref. [31] . it was shown that these phantom 
wormholes may be theoretically constructed by vanishing amounts of ANEC violating 
phantom energy, and the stability regions, in the spirit of this paper, were further 
explored in Ref. [2H]. 



5.2. Stability regions 
The surface stresses, eqs 
1 



for this particular case takes the form 



a 



V 



Ana 



1 / 1 



1 h or 



— + a 2 + ad 



6(a) 



+ a 



r,2 



:i + a&) 1 



a + ad — 



a 2 (b-b'a) 
2(a-6) 



in a 



1 _ b M + a 2 

a 



(69) 



(70) 



a v u 

The conservation identity provides us with a' = — - (a + P) + H, with 5 given by 

6'(a)a — 6(a) 



"8?r' 



6(a) 



+ 2 1 



6(a) \ $'(a) 



Equation (|3*U)l. for this particular case is given by 



a 



Ana 2 



^ + a 2 - aa 

a 

1 - + a 2 



1 _ 5^1 + m + a 2 _ ^ 



1 - ^ + a 2 

a 



which at the static solution simplifies to 

3M i 



a'{a ) 



Anal 



1 JiM 

Q.0 



1 _ 2M 



3fe(ap) _|_ b'(ao) 
2a 2 

_ b( a o) 
ao 



(71) 



(72) 



(73) 



This relationship will play an important role in determining the stability regions dictated 
by the master equation, equation (fiSjl. 

The equation of motion is of the form d 2 + V(a) = 0, and taking into account the 
factors of eqs. (f3*H|) - (|inj) . which for this case are given by 



1 



6(a)/2 + M 



and 



G 



M-b(a)/2 



the potential takes the form 
V{a) = 1 



a Aa 2 



(m 



6(a) 
2 



mi 



(74) 



(75) 



Note that the functions T and given by eqs. 



and 



respectively, are 



completely determined by the factors F, G and m s = Ana 2 a. To determine T, given by 



Stability analysis of dynamic thin shells 



13 



equation J33J), and consequently 0, given by equation ()48j). one uses the radial derivative 
of equation (|7Tj). and eqs. (|o^j) - (f7nj) . evaluated at the static solution. One may now 
model the wormhole geometry by choosing specific values for the form and redshift 
functions, and consequently determine the stability regions dictated by the inequalities 
PH|) - (|5Uj) . As the explicit form of 9 is extremely lengthy, so that as in Ref. I2ZH2HI, we 
find it more instructive to show the stability regions graphically. 

5. 3. Specific form and redshift functions 

We shall in this section consider various choices for the form function, namely b(r) = r 
and b(r) = r^/r, and deduce the respective stability regions. We shall verify that the 
latter may be significantly increased by considering appropriate choices for the redshift 
function. In this section, we shall relax the condition that the surface energy density be 
positive, as in considering traversable wormhole geometries, one is already dealing with 
exotic matter. 

5.3.1. b(r) = r . Firstly, consider the particular case of b(r) = r and $'(a) = $"(a) = 
0. The redshift function can either be constant, or have the following general form 
= J2 n >3 c n(r — a) n , so that the condition $'(a) = < ^ > "( a ) = is verified. In the 
analysis that follows below, only particular cases shall be analyzed. Note that the factor 
$"(a) of the model enters into S'. As we are matching the interior wormhole solution 
to an exterior vacuum spacetime, then there is no need to impose the condition of an 
asymptotically flat spacetime. For instance, this notion is reflected in the choice of the 
redshift function given by $(r) = (r — a) 3 /rQ. However, if one is tempted to model an 
asymptotically flat wormhole geometry, in the absence of a thin shell, then one may 
use, for instance, the choice of the redshift function given by <3>(r) = r (r — a) 3 /r 4 . In 
the analysis that follows we shall separate the cases of r < 2M and r > 2M, which 
corresponds to a positive and negative surface energy densities. 

For r < 2M, we verify that m s (a ) > and cr'( a o) < 0, and consequently a'm s < 0, 
so that the stability regions are dictated by inequality (joT?|) . The stable equilibrium 
configurations are shown below the surfaces of the plots depicted in figure El In the left 
plot, we have considered the case of $(r) = 0. Now, it is interesting to note that one 
may increase the stability regions by adequately choosing a specific redshift function. To 
illustrate this, consider for simplicity the choice of $(r) = r /r. Note that qualitatively 
the stable equilibrium regions, for this increased for high values of a/M as 

ro — > 2M, relatively to the $(r) = case. It is possible to consider more complicated 
cases, however, this shall not be done here. 

For ro > 2M, we verify that m s (ao) < and er'(ao) > 0, and the stable equilibrium 
configurations are also dictated by inequality (|5Ujl. We have considered the specific cases 
of $'(a) = $"(a) = and $(r) = —r^/r, represented in the left and right plots of figure 
H3 respectively. For both cases the stability regions are also given below the surfaces 
represented in the plots. Note that the stability regions are significantly increased, for 
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Figure 2. The stable equilibrium regions for the case of ro < 2M, with b(r) = ro, are 
represented below the surfaces. We have considered specific choices for the redshift 
function, namely, $(r) = and $(r) = ro/r, which are represented in the left and 
right plots, respectively. Note that, qualitatively, the stable equilibrium regions, for 
$(r) = ro/r, are increased relatively to the $(r) = case. See the text for details. 

the $(r) = — r /r case, as r — > 2M, as can be readily verified from figure EJ However, 
the stability regions for the $'(a) = $"(a) = are considerably greater for low values of 
M/r and of a/Vp. Thus, the message that one can extract from this analysis is that one 
may model stable wormhole geometries, by adequately choosing the redshift function. 




Figure 3. The stable equilibrium regions for the case of ro > 2M, with b(r) = ro, are 
represented below the surfaces. We have considered specific choices for the redshift 
function, namely, &'(a) = $"(a) = and <f>(r) = —ro/r, which are represented on 
the left and right plots, respectively. The stability regions are significantly increased, 
for the $(r) = —ro/r case, as ro — > 2M. However, the stability regions for the 
$'(a) = $"(a) = are greater for low values of M/ro and a/r . 

5.3.2. b(r) = r\jr. Consider the Ellis drainhole [HI], with b(r) = r\jr and $(r) = 0. 
However as in the previous case, this analysis can be extended to an arbitrary redshift 
function in which $'( a ) = < ^ ) "( a ) = 0- Recently, in an interesting paper, it was shown 
that this geometry can also be obtained with 'tachyon matter' as a source term in the 
field equations and a positive cosmological constant [35] . 
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To determine the stability regions of this solution, we shall separate the cases of 
b(a ) < 2M and 6(ao) > 2M. From equation (JBTJj) and the definition of m s = Ana^a, 
this corresponds to m s > and m s < 0, respectively. 

For b(a ) > 2M, i.e., for a negative surface energy density, and using the form 
function considered above, we need to impose the condition r > 2M, so that the 
junction radius lies outside the event horizon, a > 2M. Thus, the junction radius lies 
in the following range 

r < a < — . (76) 
o o 2M v ) 

We verify that m s a' < 0, so that according to the inequality (JHUJ) . the stability regions lie 
below the curves depicted in figure EJ We have considered the cases of r /M = 2.2 and 
ro/M = 3, where the ranges of ao are given by 2.2 < a^/M < 2.42 and 3 < a^/M < 4.5, 
respectively. Note that the values of r] are negative in the stability regions, however, 
by increasing the parameter r$/M, for $(r) =0, the values of rjo become less restricted, 
so that the stability regions are increased. 

Now, considering a fixed value of r /M and for $(r) = 0, considering the choice 
of $(r) = ro/r, we have certain regions for which the stability increases. We have 
also included the specific choice of $(r) = r/r , for instructive purposes, representing a 
specific case for which the stable equilibrium configurations decrease. 




a 



Figure 4. We have defined a = a^/M and considered the specific cases of rg/M = 2.2 
and tq/M = 3. The stability regions are depicted below the curves. For low values 
of ro/M, we verify that considering the choice of $(r) = ro/r, the stability regions 
increase relatively to $(r) = 0. For instructive purposes, we have included the case of 
<I>(r) = r/ro, for which the stability regions decrease. See the text for details. 



For b(ao) < 2M, i.e., for a positive surface energy density, and using the form 
function considered above, we shall separate the cases of r < 2M and r > 2M. 

Firstly, considering the case of r > 2M, we verify that the junction radius lies in 
the following range 



a ">m- < 77 > 
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For this specific case, a' possesses one real positive root, R, in the range of equation 
(I77jl . signalling the presence of an asymptote, a'\n = 0. We verify that a' > for 
Tq/(2M) < ao < R, and o' < for an > -R. Thus, the stability regions are given by 



Vo > @, 
Vo < 6, 



if 
if 



2M 
a > R. 



< ao < R, 



(78) 
(79) 



Consider the particular cases of r /M = 2.2, so that a /M > 2.42, and r /M = 3, 
where a /M > 4.5. The asymptotes, a'\n = 0, for these cases exist at R/M ~ 3.285 
and R/M ~ 6.50265, respectively. These cases are represented in figure El We verify 
that, for $(r) = 0, and considering a fixed high value of r /M, the specific choice of 
$(r) = ro/r increases the regions of stability. However, for low values of the parameter 
Tq/M, and for low a^/M, we need to be more careful, as certain regions for the specific 
case of $(r) = r /r, decrease the stable equilibrium configurations of the solution, as is 
qualitatively represented in the right plot of figure El 
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Figure 5. We have defined a = ao/M, and considered the cases of r /M = 3 
and ro/M = 2.2. The stability regions are depicted in the plots, (i) corresponds 
to $(r) = ro/r; (ii) corresponds to 4>(r) = 0. Note that one may increase the stability 
regions by choosing an appropriate redshift function. See the text for details. 



Considering the case of r$ < 2M, we verify that the junction radius lies in the range 
a > 2M. We have m s a' < 0, so that the stability region lies below the curves depicted 
in figure El Adopting a conservative point of view, we verify that for this case of a 
positive surface energy density, stability regions exist well within the range < 770 < 1. 
The stability regions are increased by increasing the value of ro/M, with $(r) = 0, and 
is further increased by considering the specific choice of $(r) = r /r. 



6. Conclusion 



Thin shells play an extremely fundamental role in general relativity, and have numerous 
applications ranging from the study of gravitational collapse to the Randall-Sundrum 
brane world scenario [36 . In this paper, we have analyzed the stability of generic 




Figure 6. We have defined a — ciq/M, and considered the specific choices of 
ro/M = 0.5 and r = /M = 1.25. The stability regions are depicted below the curves. 
Note that the stability regions for $(r) =0 are increased by increasing the parameter 
ro/M, and further significantly increased by choosing an appropriate redshift function, 
for instance, <f>(r) = r /r. See the text for details. 



spherically symmetric thin shells to linearized perturbations around static solutions. 
We have included the momentum flux term in the conservation identity, deduced 
from the "ADM" constraint and the Lanczos equations. Following the Ishak-Lake 
analysis, a master equation dictating the stable equilibrium configurations was deduced. 
Considering the transparency condition, we studied the stability of thin shells around 
black holes, showing that our analysis is in agreement with previous results. Applying 
the analysis to traversable wormhole geometries, considering specific choices for the form 
function, we deduced stability regions, and found that the latter may be significantly 
increased by considering appropriate choices for the redshift function. 

We have analyzed the stability by defining a parameter t], so that one does not have 
to define an equation of state of the stresses involved. Normally, y/rj is interpreted as 
the speed of sound, so that the requirement that the latter does not exceed the speed of 
light is naturally < r\ < 1. Although we have imposed this condition for the analysis 
concerning thin shells around black holes, this definition cannot be naively applied to 
stresses that violate the null energy condition, i.e., "exotic matter". Therefore, we 
have relaxed the range < r\ < 1, when considering the stability analysis for thin 
shells around traversable wormholes, although stability regions do exist for the referred 
interval. It is also interesting to note, as emphasized in Ref. |2S], that thin shells (or 
domain walls) arise in models with spontaneously broken discrete symmetries in field 
theory. Note that these models involve a set of real scalar fields (pi, with a Lagrangian of 
the form C — \ {d^i) 2 — V(<p), where the potential V(<p) has a discrete set of degenerate 
minima. Thus, by suitably choosing (pi and V((p), one may obtain the dynamically stable 
thin shells analyzed in this work. 
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